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Abstract 

We analyze recent data from the Babar and Belle Collaborations on the X(3872) resonance in 
the J /if) 7r + 7r~ and D°D°it decay channels, taking careful account of the universal features of 
an S-wave threshold resonance. Because the line shapes for such a resonance are not integrable 
functions of the energy, the resonance parameters depend on the prescriptions used to define them. 
In recent experimental analyses of the D°D°tt channel, an event near the D*°D° threshold was 
assumed to come from D*°D° or D°D*° and was therefore assigned an energy above the threshold. 
Taking this effect into account, our analysis of the D D°ir data gives a mass for the X(3872) 
that is below the D*°D° threshold. Our analyses in both the J/ijj^ir^ and D°D°ir channels 
are consistent with the identification of the X(3872) as an extremely weakly-bound charm meson 
molecule. 

PACS numbers: 12.38.-t, 12.39.St, 13.20. Gd, 14.40. Gx 
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I. INTRODUCTION 



In August 2003, the Belle Collaboration discovered a new cc meson that they named the 
X(3872) [1]. This marked the beginning of a new era in cc meson spectroscopy in which 
discoveries at the B factories have tripled the number of known cc mesons above the charm 
meson pair threshold. Of these new cc mesons, the X(3872) remains the one for which 
by far the most experimental information is available. Still the nature of this state is not 
universally recognized in the high energy physics community. 

There are two crucial pieces of experimental information that determine the nature of 
the X(3872) unambiguously. One is its mass as measured in the J/ip7t + it~ decay mode. 
By combining the most recent measurements by the Belle, Babar, and CDF Collaborations 
[2-4], its mass is determined to be Mx = 3871.55 ± 0.20 MeV. This mass is extremely close 
to the D*°D° threshold. The energy relative to the threshold is -0.25 ± 0.40 MeV. The 
central value corresponds to a bound state with binding energy Ex = 0.25 MeV. The second 
crucial piece of information is the J PC quantum numbers. Observations of decays into J/ip 7 
and ip(2S) 7 by the Belle and Babar Collaborations [5-7] imply that X is even under charge 
conjugation. The spin and parity quantum numbers have been constrained by the Belle and 
CDF Collaborations [8, 9] from studies of the angular distributions in J/^tt + tt~ decays. 
The CDF analysis is compatible only with J PC = 1 ++ and 2~ + [9]. The possibility 2~ + 
is disfavored by the observation of the decay into ip(2S) 7 [7], because it would have to 
overcome multipole suppression. The possibility 2~ + is also disfavored by the observation 
of decays into D*°D° by the Belle and Babar Collaborations [10-12], because it would have 
to overcome angular-momentum suppression associated with the tiny energy relative to the 
D*°D° threshold. We will assume from now on that the quantum numbers of the X(3872) 
are 1 ++ . 

Given that its quantum numbers are 1 ++ , the X(3872) has an S-wave coupling to the 
charm meson pairs D*°D° and D°D*°. The closeness of the mass to the D*°D° threshold 
implies that it is a resonant coupling. This state is therefore governed by the universal 
properties of S-wave threshold resonances that are predicted by nonrelativistic quantum 
mechanics [13]. We can conclude that the X(3872) is a charm meson molecule whose con- 
stituents are a superposition of D*°D° and D°D*°. Among the universa l prope rties of this 
molecule is that the root-mean-square separation of its constituents is a/ fj.Ex/2, where [i is 
the reduced mass of the D*°D°. The tiny binding energy Ex implies a large RMS separa- 
tion, with the central value 0.25 MeV corresponding to an astonishing RMS separation of 
about 6 fermis. 

Some of the confusion regarding the nature of the V(3872) has been prompted by mea- 
surements of the D D 7t and D°D 0/ -f decay modes [10-12]. In the most recent analyses by 
the Babar and Belle Collaborations, these decay modes have been analyzed as if they were 
decays into D*°D° and D°D*° [11, 12]. The resulting energy distribution must by definition 
vanish below the D*°D° threshold and it has a peak just above the threshold. Measurements 
of the position and width of this peak have been interpreted incorrectly as measurements 
of the mass and width of the X(3872). For example, in the 2008 edition of the Review 
of Particle Physics [14], the Particle Data Group determined their average for the mass of 
the X(3872) by combining four values below the threshold from J/ipn+n^ decays with two 
values above the threshold from D D°ir and _D°D°7 decays. In the PDG average mass, 
the 3.5 sigma discrepancy between these two sets of measurements was taken into account 
by increasing the error by a scale factor of 2.5. The Particle Data Group also took as their 
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average for the decay width of the X(3872) the width of the D*°D° energy distribution 
measured by the Babar Collaboration. 

Several authors have misinterpreted the measurements of the D°D°tt° decay modes as 
evidence that the X(3872) is not a bound state with mass below the D*°D° threshold 
but instead as a "virtual state" which is unbound [15-17]. The signature for a virtual 
state associated with an S-wave threshold resonance is an enhancement in the production 
of D*°D° just above the threshold together with the absence of a resonance in D°D 7i 
below the threshold. This should be contrasted with a bound state, whose signature is a 
similar enhancement above the threshold together with a resonance below the threshold. In 
misinterpreting the D°D 7r° data as evidence for a virtual state, the authors of Refs. [15-17] 
did not take into account that a bound state of D*°D° can decay into D D tt° and D°D 0/ j 
through decays of its constituent D*° or D*°. This conceptual error was pointed out in 
Ref. [18], and an analysis that takes proper account of the bound state was carried out. 

In this paper, we carry out analyses of the recent data from the Babar and Belle Collabo- 
rations on the X(3872) resonance in the J/ip7r + 7r~ and D D°ir decay channels. We begin in 
Section II by describing the line shape of an S-wave threshold resonance in a short- distance 
decay channel, such as the J/tPit + it~ decay mode of X(3872). In Section III, we use that 
line shape to analyze the most recent data from the Belle and Babar Collaborations on the 
J/i/jit + 7t~ decay channel. We proceed in Section IV to describe the line shape of X(3872) in 
the D°D 7r° decay channel, which involves the decay of a constituent. In Section V, we use 
that line shape to analyze the most recent data from the Belle and Babar Collaborations on 
the D°D°ir decay channel. We take into account the experimental procedure that identifies 
D°D°'k events near the D*°D° threshold as D*°D° or D°D*° events above the threshold. 
Our analysis of these energy distributions, which are nonzero only above the D*°D° thresh- 
old, gives a mass for the X(3872) that is below the threshold. In Section VI, we present a 
critical discussion of previous theoretical analyses of the line shapes of the X(3872). Our 
results are summarized in Section VII. 

II. LINE SHAPE IN THE J/V> tt+tt- DECAY CHANNEL 

If there is an S-wave resonance very close to the threshold for a pair of particles with 
short-range interactions, their scattering length a is large compared to the range of their 
interaction. The line shapes associated with an S-wave threshold resonance have some 
unusual features that are not ordinarily encountered in high energy physics. The line shapes 
for decay modes that involve the decay of a constituent are different from those for all other 
decay modes, because they can proceed even when the constituents have a large separation of 
order a. For all other decay modes, the constituents must approach to within a much smaller 
distance comparable to the range of the interaction. We will refer to these two classes of 
decay modes as constituent decay modes and short- distance decay modes, respectively. In the 
case of the X(3872), the constituent decay modes are D°D 7r° and D°D 0r j and an example 
of a short-distance decay mode is J/ip7i + n~. In this section, we summarize the essential 
aspects of the line shape for short- distance decay modes. 

The line shape for a resonance near a scattering threshold is proportional to \f(E)\ 2 , 
where f{E) is the analytic continuation of the scattering amplitude in the total energy E of 
the particles in their center-of-mass frame. The universal scattering amplitude for an S-wave 
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threshold resonance has the form 



f(E) = , (1) 

-7 + V-ME + ie) 

where E is the energy relative to the threshold, \i is the reduced mass, and 7 = 1/a is 
the inverse scattering length. In the case of the X(3872) resonance, the relevant scattering 
amplitude is for D*°D° mesons in the 1 ++ channel. For quantitative applications, the 
scattering amplitude in Eq. (1) must be modified to take into account the nonzero width 
of the D*° and the existence of inelastic scattering channels for the charm mesons [18]. By 
analytically continuing the parameters in Eq. (1) to complex values, we obtain 

(2) 



f(£J) — 

-(7re + Him) + y/-2fJ,(E + iT* /2) ' 



where \i = 966.6 MeV is the reduced mass of the D*° and D°, r* = 65.5 ± 15.4 keV is the 
total width of the D*°, and Tre + Him is the complex inverse scattering length. The effects of 
the decays of D*° into D°7r° and _D°t are taken into account through r*o- The effects of the 
inelastic scattering channels for D*°D°, such as J/tpn + n~, are taken into account through 
7i m , which must be positive. The scattering amplitude in Eq. (2) should be accurate as long 
as the energy is within about an MeV of the threshold. 

For a short-distance decay channel, the only dependence of the line shape on the energy E 
is from the resonance factor \f(E)\ 2 . If Tre is positive, the line shape \f(E)\ 2 has a resonance 
peak below the D*°D° threshold. Defining the binding energy and decay width for this 
resonance is problematic, because the line shape is not that of a conventional Breit-Wigner 
resonance. Our prescriptions for the binding energy E x and the width T x are that the pole 
of the amplitude f(E) in the complex energy E is at —E x — iT x /2: 
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7-1 /re /im /o \ 

E X = —j—, (3a) 

T X = r, + . (3b) 

In the case T x <C 2E X , the shape of the resonance is approximately that of a nonrelativistic 
Breit-Wigner resonance in the region \E + E X \ <C E x . Its peak is at —E x and its full width 
at half maximum is T x , justifying the interpretation of E x and T x as the binding energy 
and decay width of the resonance. If T X /(2E X ) is not small, the variables E x and T x 
defined by Eqs. (3) have no precise physical interpretations. 

If Tre is negative, the line shape \ f(E)\ 2 has a peak very near the D*°D° threshold. In 
the limit r^ — > 0, the peak is a cusp with a discontinuity in the slope that arises from the 
square root in Eq. (2). The effect of the D*° width is to smooth out the cusp. In this case, 
the variables E x and T x defined by Eqs. (3) specify the location of a pole on the second 
sheet of the complex energy E. Thus they have no simple physical interpretations. 

The binding energy E x and the width T x can not be measured directly, because they are 
defined in terms of the analytic continuation of the scattering amplitude f(E) to complex 
values of the energy E. An alternative pair of variables that can in principle be measured 
directly are the position E max of the peak in the line shape and its full width at half-maximum 
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Tfwhm- The position E mSLX of the peak satisfies 

\ 1/2 



+ 7im (/^max + ry4 + /^max) 7 = . (4) 

The full width of the line shape at half- maximum is given by rf w h m = E + — E_, where E± 
are the two solutions of 

l/(£±)| 2 = ^(£max)| 2 . (5) 

If 7 re > 0, the solutions for E max , E + , and E_ can be expanded in powers of r* - The 
expansions for E max and rf w hm are 

£ max = ~ (il + ^r* ) + ll ' 3 4 lL (^ ) 2 + ■■) , (6a) 

^f 1 \ Ire 47 re J 

r fwhm = (4 7re7im + 2(//r, ) + ^f^f^ ^T^) 2 + ..Y (6b) 

A* V TrevTre Tim/ / 

The normalization of the line shape of X in a short-distance decay mode F produced 
by the decay of B + into K + + X is proportional to the product of the branching fractions 
for B + — > X + + X and X — > F. It is convenient to introduce a compact notation for the 
product of these two branching fractions: 

(BB) F = B[B + -> X+ + X] £[X -> F] . (7) 

Defining these branching fractions is problematic, because the line shape for an S-wave 
threshold resonance is not an integrable function. Since \f(E)\ 2 decreases as 1/\E\ for large 
\E\, the integral of \f(E)\ 2 over E depends logarithmically on the endpoints. This implies 
that this product of branching fractions cannot be defined uniquely in terms of an integral 
over the line shape. The numerical value of {BE) p will inevitably depend on the prescription 
used to define it. Our prescription is that the normalized line shape for B ± — > K ± + F is 

^ r|B+1(Be ) F ^, (8) 



where the energy-dependent factor is 



dv SD /x 2 r 



X 



dE 27r( 7r 2 e + Tfm) 



\f(E)\ 2 ■ (9) 



In the case Tx 2Ex, this line shape in the region \E + Ex\ <C Ex is well approximated 
by a Breit-Wigner resonance. The integral of dT S D/dE over this region is approximately 1, 
justifying the interpretation of (BE) f as the product of the branching fractions for B + — > 
K + + X and X -> F. If T X /{2E X ) is not small, the constant {BB) F defined by Eq. (8) 
has no precise physical interpretation. It is simply a convenient variable for specifying the 
normalization of the line shape. An alternative prescription for (BB)p could be obtained 
by integrating both sides of Eq. (8) over a chosen interval of the energy E in the threshold 
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FIG. 1: Invariant mass distribution for the J/ipTr + ir~ decay channel measured by the Babar 
Collaboration [2]. The data are the number of events per 5 MeV bin. The inverse scattering 
lengths 7 re + i^i m for the two fits are 38.8 MeV (dashed line) and (13.6 + 15. 5i) MeV (solid line). 
The vertical line is the assumed D*°D° threshold at 3871.8 MeV. 



region, such as —2E X to 0. However the numerical value of (BB)p would depend on the 
choice of the endpoints of the interval. 

Although the product of branching fractions depends on the prescription, the ratio of 
{BB)f for two short-distance decay modes F is independent of the prescription. Choosing 
one of the final states to be J/ip7r + n~ and using Eq. (7), the ratio is 

(BB) F = B\X -> F] 
(BB) J/il7r+7r - B[X^J/^n+n-] ■ 1 ] 

The ratio on the right side of Eq. (10) is the conventional branching ratio for decays of X 
into those states. This ratio is well-defined for any short- distance decay mode F, despite 
the fact that a prescription is required to define the products of branching fractions on the 
left side of Eq. (10). 



III. ANALYSIS OF THE J/rjj tt+tt- DECAY CHANNEL 

In this section, we analyze recent data from the Belle and Babar Collaborations on the 
line shape of the X(3872) in the J/ ip 7r + 7r~ decay mode [2, 3]. We consider the invariant mass 
distribution for J/if)7r + 7i~ in the interval from 3820 MeV to 3920 MeV. For our two data 
samples, the total number N B g of B + B~ and B°B° events accumulated and the number of 
candidate events for the decay of B ± into K ± + J/ipTt + T{~ are as follows: 

• Babar Collaboration [2]: N bS = 4.55 x 10 s , 471 events in 20 bins of width 5 MeV, 
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FIG. 2: Invariant mass distribution for the J/i^tt^tt decay channel measured by the Belle Col- 
laboration [3]. The data are the number of events per 2.5 MeV bin. The inverse scattering lengths 
7re + «7im for the two fits are 47.5 MeV (dashed line) and (38.4 + 12.0i) MeV (solid line). The 
vertical line is the assumed D*°D° threshold at 3871.8 MeV. 



• Belle Collaboration [3]: N BB = 6.57 x 10 8 , 606 events in 40 bins of width 2.5 MeV. 

The data are shown in Figs. 1 and 2. The vertical error bar in a bin with n events is ^/n. 
The horizontal error bar indicates the width of the bin. 

We take the theoretical line shape for the energy E of J/ipn + ir~ relative to the D*°D° 
threshold to be given by Eqs. (8) and (9). This line shape will be accurate within about an 
MeV of the threshold. We assume that the dominant contributions to the signal come from 
this threshold region. If this is the case, then a line shape that remains accurate over a larger 
energy interval would give a better approximation only to contributions that are negligible. 
To obtain a line shape that remains accurate within about 10 MeV of the threshold, it is 
necessary to take into account the effects of the charged charm meson pairs D* + D~ and 
D + D*~ , as discussed in Section VI. 

To predict the number of events in a given bin of invariant mass, we need to take into 
account the background and the energy resolution of the experiment. The resolution must 
be taken into account because the line shape varies dramatically over an energy scale smaller 
than the energy resolution. The predicted number of J/ip7T + 7i^ events in an energy bin of 
width A centered at ~E% can be expressed as 



N; = 2N 



i — ^BB 



pEi+A/2 poo rp 

(BB)j/^- dE' d ER(E',E)—^ + C hg A 

JEi-A/2 J-oo atj 



where C*b g is the background under the line shape dTsD/dE. Our invariant mass interval 
3820-3920 MeV is narrow enough that we take the background term C*b g to be a constant 
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TABLE I: Results of our analyses of the data for — ► + J/t/j ir + Ti~ . The four rows correspond 
to analyses using either the Babar data [2] or the Belle data [3] and either setting y im = or using 
7i m as a fitting parameter. All entries are in units of MeV, except for (BB) J /^ n + 7r -, which is in 
units of 10~ 6 . 



independent of E. The experimental resolution is taken into account through the convolution 
with the Gaussian resolution function: 

R(E>, = exp(-(£' - £) 2 /(2<x 2 )) • (12) 

\/ Itxo 

We follow Ref. [17] in taking the width of the Gaussian to be the same energy-independent 
constant for both experiments: a = 3 MeV. 

We assume that the number of events in each bin of the smeared J/ipn + 7i~ energy E' 
has a Poisson distribution whose mean value is given by JVj in Eq. (11). We fix the D*°D° 
threshold at 3871.8 MeV and the D*° width r„, at 65.5 keV. The fitting parameters are 7 re , 
7i m , (BB)j/^ 7T + n -, and C^ g . We determine the best fit to these parameters by maximizing 
the likelihood for the observed distribution. For both the Belle and Babar data sets, we 
carry out two fits, one with 7 im = and one with 7 im as a fitting parameter. The results of 
our four analyses are presented in Table I. The error bars on 7 re and 7i m are determined by 
varying these parameters while keeping (BB)j/^ n + n - and Cb g fixed at their central values. 
For 7 im = 0, the error bars on 7 re give the interval within which log(Likelihood) differs 
from its maximum value by less than 1/2. If 7 im is treated as a fitting parameter, the 
error bars for 7 rc and 7 im specify the smallest rectangle that contains the error ellipse in 
which log(Likelihood) differs from its maximum value by less than 1/2. The error bars on 
(BB)j/ t (, 7T + 7T - are determined by varying this parameter and Cb g while keeping 7 re and 7i m 
fixed at their central values. 

In Table I, we also give the calculated values of the position of the resonance and its 
width using two different prescriptions for the parameters. The values of —E x and Y x were 
calculated using Eqs. (3). The values of E max and r fwhm were obtained by solving Eqs. (4) 
and (5). The uncertainty of ±0.36 MeV in the energy of the D*°D° threshold is taken into 
account as an additional statistical error in —Ex and in E max . The uncertainty of ±15.4 keV 
in r* is taken into account as an additional statistical error in Tx and in Ff w h m - In Table I, 
there are significant differences between the values of —Ex and E max for the fits in which 7 im 
is used as a fitting parameter. All four fits give values of (BB)j/ 1 p 7T + 7T - that are consistent 
to within the errors and approximately equal to 10~ 5 . 

The fits to the Belle data give parameters 7 rc and 7 im with smaller error bars than the fits 
to the Babar data. In Figs. 3 and 4, the unsmeared line shapes of X(3872) in the J/^7t + tt^ 
decay channel corresponding to the central values of the two fits to the Belle data are shown 
as solid lines. The line shape in Fig. 4 from using 7i m as a fitting parameter is wider than 
that in Figs. 3 from setting 7i m = 0. Both line shapes are much narrower than the smeared 
line shapes shown in Fig. 2. Thus most of the observed width can be accounted for by the 
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FIG. 3: Line shapes of X(3872) for j ve + i^i m = 47.5 MeV. The curves are the line shape in 
J/ipTT + iT~ (solid line), the line shape in D°D°'K (dashed line), and the D*°D° energy distribution 
(dash-dotted line). The two line shapes have been normalized so the resonances below the threshold 
have the same peak height. 



experimental resolution. 



IV. ENERGY DISTRIBUTIONS FOR THE D°D°tv° DECAY CHANNEL 

In this section, we summarize the essential aspects of the line shape of the X(3872) 
in the D°D°7i° channel. We also determine the energy distribution that follows from the 
identification of D D°ir events with energy near the D*°D° threshold with D*°D° and 
D°D*° events above the threshold. 

In the decay B + — ► K + + D D tt°, the momentum distributions for D D tt° near the 
X(3872) resonance can be calculated from the sum of the two diagrams in Fig. 5. The open 
dot represents the B + — > K + transition which creates a D*°D° or D°D*° at short distances. 
The double line represents the exact propagator for the resonant superposition of D*°D° 
and D°D*°, whose dependence on the total energy E of D°D 7r° is given by the scattering 
amplitude f(E) in Eq. (2). In the propagators for the virtual D*° and D*°, the width r* 
must be taken into account. The coupling of the ir° to the charm mesons is linear in the 
pion momentum. The differential distribution in the total energy E and in the momenta 
Pd, Pd-i an d Ptt of the D°, D°, and 7i~° has the form 



dY^\f{E)\ 2 V 2 



+ 



1 



p 2 D - 2/iE - z/ir* p 2 D - 2fiE - ifiT^Q 



d<S> D£)7T dE . (13) 



The differential 3-body phase space d&pjj^ includes a delta function that relates the energy 
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FIG. 4: Line shapes of X(3872) for 7re + i^ im = (38.4 + 12.0*) MeV. The curves are the line 
shape in J/t/jir + TT~ (solid line), the line shape in D D°ir (dashed line), and the D*°D° energy 
distribution (dash-dotted line). The two line shapes have been normalized so the resonances below 
the threshold have the same peak height. 





FIG. 5: Diagrams for the production of D°l) 7r . The open dot represents the B — > K transition 
that creates D*°D° or D°D*° at a short-distance scale. The double line represents the propagation 
of the resonant linear combination of the pair of charm mesons. The two diagrams involve either 
a virtual D*° (left diagram) or a virtual D*° (right diagram). 



E and the three momenta: 



E = —Sn*Dn + 



Pi 



+ 



2M D o 2M D o 2m^ 
where 5d*Dtt is the energy released in the decay of D*° to D°7r°: 

5 D * DlT = M D ,o - M D o - m^o = 7.14 ± 0.07 MeV. 



(14) 



(15) 
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The 3-body phase space can be reduced to a 2-dimensional integral over pn and pq\ 



2m 7r o 

d$DDn = JT^PDdpDPDdpD- (16) 

The integration region is 

— ~ < 0D*Dtt + < —Z 1 , [i-l) 

2fi D7T m w o 2fi D7T m^o 

where /inn = M D om n o / (M D o + m n o) is the reduced mass of D° and n°. 

Voloshin has used the diagrams in Fig. 5 to predict the momentum distributions for 
D°D°ir in decays of the X(3872) resonance [19]. His result is proportional to the right side 
of Eq. (13) with the resonance factor \ f(E)\ 2 omitted, r*o set to 0, and the energy E in 
the propagators replaced by —Ex, where Ex is the binding energy of the X(3872). This is 
the appropriate momentum distribution only if the energy E is fixed at a value close to the 
peak of the resonance. In the case of a low-energy antiproton beam incident on a nucleon 
target, it may be possible to tune the center-of-mass energy to the peak of the resonance. 
However in the case of B meson decays, the X(3872) resonance is produced with a variable 
energy E. Since the experimental resolution in E is larger than the width of the resonance, 
it is necessary to take the resonance factor \f(E) | 2 into account. 

If \E\ is small compared to 2(m 7r o /M D o)8 D * D7T = 1.04 MeV, the phase space integral in 
Eq. (13) can be evaluated analytically. 1 In this case, the interference term between the 
two propagators in Eq. (13) can be neglected. The 3-body phase space integral in Eq. (13) 
reduces to 

2 1 /9„E LAI. \V2 / , xl/2 



- - ' l- f '^^>^\ Ue 2 + Y 2 JA + e\ .(18) 



\p 2 D - 2fiE - ifiT^l 2 "' DD7T 7r 2 r* V A* 
The resulting line shape has the form 



^K\f(E)\ 2 ^E 2 + Tl /4 + E y 2 . (19) 

This simple expression for the line shape was first derived by Braaten and Lu [18]. If 7 rc 
is positive, the line shape consists of a resonance associated with the bound state below 
the D*°D° threshold and a threshold enhancement above the threshold. If 7 rc is negative, 
there is a threshold enhancement above the D*°D° threshold but no resonance below the 
threshold. The position £^ n of the maximum in the line shape satisfies 

2/V £ - + r 2 /4 - 4/i£ max + 7r 2 c + 7i 2 m - 4 7 rc (» yj + Y%/A - ptE m Jj ' =0. (20) 
The solution up to corrections that are second order in T* is 

F D°BW _ ( 2 7rc + V7rc ~ Hn\ 

~ 2/i ' 1 



1 It is also necessary for to be small compared to (m n o /M d o)5d*Dk, but this is satisfied if r* ~ 

65.5 keV. 
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The normalization factor for the line shape in Eq. (19) involves the product (BB) D o£)o w o of 
the branching fractions for B + — > K + X and X — > D°D°7r° defined by Eq. (7). Defining these 
branching fractions is problematic, because the line shape in Eq. (19) is not integrable. Since 
\f(E)\ 2 decreases as 1/\E\ for large \E\, the integral of the line shape in Eq. (19) increases as 
the square root of the upper endpoint. This implies that the product of branching fractions 
cannot be defined uniquely in terms of an integral over the line shape. The numerical value 
of (BB) D oQo n o depends inevitably on the prescription used to define it. Our prescription is 
that the normalized line shape for B^ — > K ± + D°D°ir° is 

1/2 



dv -t\r+\(bb\ - dFsD ( _V£+£kt±^- I 
- = r [B] pw,-^^=_| , (22) 

where dtso/dE is the short- distance line shape in Eq. (9) and Ex is the binding energy given 
by Eq. (3a). The last factor in Eq. (22) reduces to 1 at E — —Ex- In the case Tx <C 2Ex, 
the line shape in the region \E + Ex\ <C Ex is approximately a Breit-Wigner resonance. 
The integral of the right side of Eq. (22) over this region is approximately r[£> + ] (£>£>) d°d°tt j 
justifying the interpretation of (BB)f as the product of the branching fraction for B + — > 
K + + X and X -> D°D n°. If Tx/(2E X ) is not small, the constant (BB) D o D w defined 
by Eq. (22) has no precise physical interpretation. It is simply a convenient variable for 
specifying the normalization of the line shape. 

To compare with the energy distribution for D*°D° measured in the B factory exper- 
iments, we must take into account how these energy distributions are measured. Three 
particles identified as D°, D°, and n can be considered as candidates for either a D*°D° 
event or a D°D*° event. If the measured invariant mass of D°ir° is close enough to the mass 
of D*° (within 10 MeV for Babar [11], within 6 MeV for Belle [12]), it is identified as a D*°. 
The constraint that the invariant mass of D°tt° is equal to Md*o is then used to sharpen 
the resolution of the measured momenta. If the D° and n° are produced by the decay of 
a constituent D*° from the bound state X(3872), their invariant mass will be smaller than 
M D *o by approximately the binding energy E X - This information about the binding energy 
is discarded when the D V is constrained to come from the decay of a D*°. If the momenta 
of the D°, D°, and 7r° in the D°D°7r° rest frame are pn, Pd, an d Pw, their total energy E 
relative to the D*°D° threshold is given in Eq. (14). If the D 7i° is identified as a D*° in 
the experimental analysis, the inferred energy E cxp of the D°D°ir° relative to the D*°D° 
threshold is 

(pp+Pn) 2 Pp = Pi m) 
2M D *o 2M D o 2/1 ' 1 1 

Similarly, if the D°7r° is identified as a D*° in the experimental analysis, the inferred energy 
E exp of the D°D°n relative to the D*°D° threshold is p 2 D /(2fi). If the D°ji° and D°n both 
have invariant mass close enough to M D *o to be identified as D*° and D*°, the one whose 
invariant mass is closest is constrained to be a D*° or D*°. Thus the inferred energy E exp 
of a D°D 7v° event that is identified as either D*°D° or D°D*° is 

E cxp - E < (24a) 

- max (p 2 D>Pp) E pI + Pd mh) 

2/i 4/x 1 1 



12 



We will refer to E cxp as the D*°D° energy. It is the distribution in this variable that should 
be compared with the experimental energy distribution. 

If \E\ is small compared to 2(m^-o /M D o) S d*d-k = 1-04 MeV, the line shape in the D°D°ir 
channel is given by the analytic expression in Eq. (22). In this case it is also possible to obtain 
an analytic expression for the distribution in the variable E cxp . For \E\ <C 2(m 7r o / 'M d o)5d* o-k, 
the larger of the momenta po and p^ is approximately (2fi Dw 5D* d-k) 1 ^ 2 , s° E cxp is given by 
Eq. (24a). The distribution in E exp can then be expressed as 



dr _ r[B+}(BB) DODOnO r* . /2 r°° 



< 2 / dE^j= =^™, (25) 



dE cxp V2n{^Ej+T%/4-E x ) 1/2 ° XP J-oo dE \E cxp — E — iT* /2\ 2 

where Tso/dE is given in Eq. (9) and Ex is given in Eq. (3a). The normalization is consistent 
with that in Eq. (22), as can be verified by integrating over E exp using the integration formula 

The integral over E in Eq. (25) can be evaluated analytically by deforming the integration 
contour into the upper half-plane and picking up the contributions from the two poles and 
the branch cut. The resulting expression for the integral over E in Eq. (25) reduces to 



j — < 



, „ dT SD 1 n 2 T x 

dE 



dE |£ cxp -£-TV2| 2 27rr, |7| 2 

2i 1 2 Q F(-i 1 2 ) 



( 7 2 _ % 2 + 2i7 2 )(7 2 + 2/iE exp )(7 2 + 2fiE exp + 2* 7o 2 ) 

22 7o 2 [F{-i^ 2 - 2 7o 2 ) - 27r( % + lyT^) (7 + V7 2 = 2^ 2 )] 
( 7 2 _ 7 2 + 2i7 2) (7 2 + 2/i£ cxp )( 7 2 + 2^ cxp - 2? 7o 2 ) 

F{2iyE exp - 2 7 2 ) - 27r( 7 , + i^/2fiE cxp ) (7 - i^/2fiE exp + 2i 7o 2 ) 
( 7 2 + 2/i£ cxp )( 7 2 + 2/^ exp + 2^ 7o 2 ) 
F(2tfiE cxp ) 



( 7 2 + 2/i£ cxp )( 7 2 + 2fiE cxp - 2* 7o 2 ), 
where 70 = (fiT^) 1 / 2 , 7 = 7 re + ^7i m , 7* = 7™ - ?7im, and the function F(z) is 

= « J -2(2 + 2 7 2 ) 2tt 7 - AV^iz log - — — ^ i 'OJ. 

v \ 270 



(27) 



(28) 



This function has a square-root branch point at z — — 2 / -/q, but despite the factors of \fz it 
has no branch point at z — 0. Although it is not manifest, the expression on the right side 
of Eq. (27) is real-valued. 

In Figs. 3 and 4, the solid lines are the line shapes in the J/ipir + i[~ decay channel for 
7rc + ^7im = 47.5 MeV and (38.4 + 12. Oi) MeV, respectively. For comparison, the line shapes 
in the D°D°7i° decay channel and the D*°D° energy distributions are also shown as dashed 
and dash-dotted lines, respectively. In each figure, the curves are normalized so that the 
resonances below the threshold have the same maximum values. In both figures, the D D n 
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FIG. 6: Energy distribution for D*°D° and D°D*° measured by the Babar Collaboration [11]. The 
data are the number of events per 2 MeV bin. The inverse scattering length 7 re + «7i m for the fit 
is 67.7 MeV. 



line shape has a peak below the D*°D° threshold corresponding to the X(3872) resonance 
and a second peak above the threshold corresponding to a threshold enhancement in the 
production of D*°D° and D°D*°. The position and width of the resonance peak is close to 
that for the J/ip7t + iT~ line shape. The D*°D° energy distribution, which vanishes below 
the threshold, has a peak above the threshold whose width is considerably larger than the 
width of the resonance. Thus a measurement of the position and width of the peak in the 
D*°D° invariant mass distribution should not be interpreted as a measurement of the mass 
and width of the X(3872). 



V. ANALYSIS OF THE D°D 7v° DECAY CHANNEL 

In this section, we analyze recent data from the Belle and Babar Collaborations on the 
line shape of the X(3872) in the D°D°7[° decay mode [11, 12]. We consider the energy 
distribution for D*°D° and D°D*° in the interval from to 100 MeV. For our two data 
samples, the total number N BB of B + B~ and B°B° events accumulated and the number of 
candidate events for the decay of B ± into + D* D°(D°D* ) £11*6 cLS follows: 

• Babar Collaboration [11]: N BB = 3.83 x 10 8 , 172 events in 50 bins of width 2 MeV, 

• Belle Collaboration [12]: N BB = 6.57 x 10 s , 171 events in 50 bins of width 2 MeV. 

The data are shown in Figs. 6 and 7. The vertical error bar in a bin with n events is \fn. 
The horizontal error bar indicates the width of the bin. 
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FIG. 7: Energy distribution for D*°D° and D°D*° measured by the Belle Collaboration [12]. The 
data are the number of events per 2 MeV bin. The inverse scattering length ~f re + ijim for the fit 
is 9.99 MeV. 



We take the distribution in the D*°D° energy E cxp to be given by Eqs. (25), (27), and 
(28). This energy distribution should be accurate within about an MeV of the threshold. We 
assume that the dominant contributions to the signal come from this threshold region. To 
predict the number of events in a given bin of invariant mass, we need to take into account 
the background and the energy resolution of the experiment. The predicted number of 
D*°D° events in an energy bin of width A centered at E{ can be expressed as 



rEi+A/2 poo 

2N BS dE' dE exp R(E',E exp )El£ 

JEi-A/2 JO 



-A/2 



V2n(^/E' x + T'i /A-E x ) 



1/2 



dE 



dY 



SD 



dE \E cxp -E -2lV2| 2 



(29) 



where C^ g takes into account the background. Our energy interval 0-100 MeV is narrow 
enough that the background contribution to the distribution in E exp can be taken as a 

constant Cb g multiplied by Ellp, which is the energy dependence of the D*°D° phase space. 
The experimental resolution is taken into account through a convolution with a Gaussian 
resolution function with an energy-dependent width a(E exp ): 



R(E', E, 



1 



exp J 



2tt(j(E i 



exp(-(£' 



£exp)7(M£exp) 2 )) 



(30) 
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0.066 ±0.015 
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U.UDD_ 0015 


9 07+O.6I 
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66+ 025 
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Belle 
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q qq+1-99 
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n +0.98 
u -0 


n 0029 +a0007 
u.uu^»_ 0006 


-n 052+°- 014 

U.UOZ_ Q 023 

-n 052+°- 014 

u.uoz_ 038 


0.066 ±0.015 

066+ 025 
u.uoo_ 0015 


-0 056+ 013 

U.UOD_ 022 
U.UOO_ 0037 


0.066 ± 0.015 
66+ 026 

U.UOD_q 015 



TABLE II: Results of our analyses of the data for — > if 1 * 1 ±.D Z) 7r . The four rows correspond 
to analyses using either the Belle data [12] or the Babar data [11] and either setting 7; m = or 
using 7i m as a fitting parameter. All entries are in units of MeV, except for (BB) D oQo n o, which is 
in units of 10~ 6 . 



We follow Ref. [17] in taking the width for both experiments to be the same energy- dependent 
function: 

a(E cxp ) = ^(0.031 MeV)£ cxp . (31) 

This may be too crude a model for the effects of the experimental resolution in this channel, 
but we will use it for illustrative purposes anyway. 

We assume that the number of events in each bin of the smeared energy E' has a Poisson 
distribution whose mean value is given by Ni in Eq. (29). We fix the D*° width r*o at 
65.5 keV. The adjustable parameters are 7 re , 7i m , (£>£>) d d 7t > an d Cbg- We determine the 
best fit to these parameters by maximizing the likelihood for the observed distribution. For 
both the Belle and Babar data sets, we carry out two fits: one with 7 im = and one with 
7 im as a fitting parameter. The results of our four analyses are presented in Table II. The 
error bars are determined in the same way as those in Table I, except that the uncertainty of 
±0.36 MeV in the D*°D° threshold energy does not enter because the experimental energies 
were measured relative to this threshold. 

In the fits to the Babar and Belle data with 7 im treated as a fitting parameter, the 
maximum likelihood is obtained for 7 im = 0, which is the smallest possible physical value. 
This suggests that our model for the experimental resolution in Eq. (31) may provide too 
much smearing of the energy distribution. The best fit to the Babar data gives a line shape 
whose peak is below the D*°D° threshold by about 2.4 MeV, which is incompatible with 
the assumption \E\ <C 1 MeV that we used to derive analytic expressions for the line shape 
and the energy distribution. The best fit to the Belle data gives a line shape whose peak 
is below the D*°D° threshold by only about 0.05 MeV, which is comparable to the width 
Fx ~ 0.07 MeV. The value of (BB)^^^ from the fit to the Belle data is about an order 
of magnitude smaller than that from the fit to the Babar data. Since the Belle fit does 
not satisfy T x 2E X , the value of (BB) d°dW should not be interpreted literally as the 
product of branching fractions. It is simply a parameter used to specify the normalization of 
the line shape in Eq. (22). The Babar fit does satisfy Tx <S 2Ex, so the value of (i3i3) D o£ ) o 7r o 
can be interpreted as the product of branching fractions. Dividing by the value 8.7 x 10 -6 
for (BB) J / 1 p 7T + 7T - from Table I, we obtain a branching ratio for D°D°-k to J/ipir + iT~ of 
approximately 0.004. This implies that short-distance decay modes account for most of the 
width Tx of the X(3872) resonance. However the value 7i m = for the best fit implies that 
the contribution to the width Tx from short-distance decay modes is negligible. A possible 
explanation for this inconsistency is that the simple model for the D*°D° energy resolution 
given in Eqs. (30) and (31) is inadequate. 
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FIG. 8: Line shapes of X(3872) for j Te + i^f\ m = 9.99 MeV. The curves are the line shape in 
J/i/jit + it~ (solid line), the line shape in D D°tt (dashed line), and the D*°D° energy distribution 
(dash-dotted line). The two line shapes have been normalized so the resonances below the threshold 
have the same peak height. 

In Fig. 8, we show the line shapes corresponding to the best fit to the Belle data for 
D°D 7i°. The line shape in the J/ip ir + 7r~ decay channel, the line shape in the D D°ir 
decay channel, and the D*°D° energy distribution are shown as solid, dashed, and dash- 
dotted lines, respectively. The curves are normalized so that the resonances below the 
threshold have the same peak height. The D°D°tt° line shape has a single peak below the 
D*°D° threshold whose position and width are close to those for the peak in the <//■?/> 7r + 7r~ 
line shape. It is this peak that should be identified with the X(3872) resonance. The 
D*°D° energy distribution, which vanishes below the D*°D° threshold, has a peak above 
the threshold whose width is considerably larger than that of the X(3872) resonance. It is 
also much narrower than the smeared energy distribution shown in Fig. 7. 

VI. CRITIQUE OF PREVIOUS ANALYSES 

In this section, we discuss how the analysis presented in this paper could be improved. 
We also point out errors and misconceptions in previous theoretical analyses of the line 
shapes of the X(3872). 

The most limiting aspect of our analysis was the use of the analytic expression in Eq. (22) 
for the line shape in the D°D 7i° decay channel. The derivation of this expression involved 
the assumption \E\ <C 2(771^0 /M D o) 5 d*d-k = 1-04 MeV. It requires most of the D°D°tt 
events to be produced with energy within an MeV of the D*°D° threshold. In particular, it 
requires the width of the X(3872) resonance to be much less than an MeV. This limitation 
can be relaxed by replacing the invariant mass distribution in Eq. (22) by the differential 
decay distribution in Eq. (13). Given our prescription for (BB) D ojjo^o in Eq. (22), the 
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normalization of the differential decay rate is determined: 
dT Y[B + }{BB) DOD ^ou 2 Y^ ( /I 



dE ~ (^i + ry4 + E x y/* \W D Jl* D « 



1/2 



dV 2 



SD 2 

X ^E P * 



+ 



p 2 D - 2/iE - ifiT^o p 2 D - 2/iE - i[ir* 



d$ DDn , (32) 



where d^ D Q n is given in Eq. (16). The integral over the momenta po and pfj must be 
evaluated numerically. 

Another limiting aspect of our analysis was ignoring the effects of the charged charm 
meson pairs D* + D~ and D + D*~. They can produce significant interference effects for 
\E\ as small as 1/4 of the 8.1 MeV splitting between the D* + D~ and D*°D° thresholds 
[20]. The effects of charged charm meson pairs were first considered by Voloshin [16], but 
there were conceptual errors in his analysis. A correct analysis was presented by Braaten 
and Lu [20]. It involves the 2x2 matrix of S-wave C = + scattering amplitudes fij(E) 
between the neutral channel D* D°/D°D* labelled by subscript and the charged channel 
D* + D~ /D + D*~ labelled by subscript 1. The parameters in these scattering amplitudes are 
the inverse scattering lengths 70 and 71 for charm mesons in the channels with isospin and 
1 in the isospin symmetry limit. The coupled-channel expressions for the D*D scattering 
amplitudes are 

-( 7o + 7l ) + 2 Kl (£) 
Joo{E) - D\E) ' ^ ' 

foi(E) = f 10 (E) = (33b) 
ME) = -^ + ^ E ] ME \ (330 

where the denominator is 

D(E) = 2 7o7l - ( 7o + 7 i)[«(£) + + 2 Ki (E)k(E). (34) 

The functions k(E) and k±(E) are 

k{E) = y/-2ij[E + iT, /2}, (35a) 
= y/^fi^E-u + iT^], (35b) 

where ^ = 968.7 MeV is the reduced mass of D* + and D~ and v = 8.1 MeV is the splitting 
between the D* + D~ and D*°D° thresholds. In Eq. (35a), r* is the energy-dependent width 
of the D*°, which has its physical value 65.5 keV at E — 0. In Eq. (35b), r*i is the energy- 
dependent width of the D* + , which has its physical value 96 keV at E — 8.1 MeV and 
decreases to 1.5 keV at E — 2.1 MeV, which is the D + D~-k° threshold. Near the D*°D° 
threshold, the scattering amplitude in Eq. (33a) reduces to the universal expression in Eq. (2) 
with 

7,. + .7 im = 27 f (36) 

(7o + 7i) - 2/tiO) 

Voloshin's first conceptual error in Ref. [16] was assuming that a B — > K transition that 
produces the X(3872) resonance must create the charm mesons in the neutral channel and 
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not in the charged channel. This implies that the amplitudes for the resonant production of 
an isospin-0 final state such as J/ifj7r + 7i^7i° and an isospin-1 final state such as J/ipit+n^ 
are proportional to f 00(E) — f 01(E) and f 00(E) + f i(E), respectively. However, since there 
is resonant scattering between the neutral and charged channels, the X(3872) resonance can 
also be produced by a B — > K transition that creates the charm mesons in the charged 
channel. Thus the amplitudes for isospin-0 and isospin-1 final states can also have terms 
proportional to f 10(E) — fn(E) and f 10(E) + fn(E), respectively. Voloshin's second con- 
ceptual error was ignoring the constraints of isospin symmetry on the amplitudes for the 
creation of charm meson pairs by the B + — > K + and B° — * K° transitions. He concluded 
incorrectly that the rates are proportional to \f 00(E) ±foi(E)\ 2 with the same multiplicative 
constant for B + and B° decays. The isospin symmetry constraints derived in Ref. [20] imply 
that the line shapes are different for B + and B° decays and they are determined by three 
independent multiplicative constants. 

Another limiting aspect of our analysis was ignoring the effects of the P-wave charmonium 
state Xci(2-P). This state has the same quantum numbers 1 ++ as the X(3872), so it also 
has an S-wave coupling to charm meson pairs D*D and DD* . Its spin symmetry partner 
X c2 (2P), which was discovered in 2006 by the Belle Collaboration [21], has a mass of about 
3930 MeV. Potential models predict the mass of the xa(2P) to be lower by 20 to 50 MeV 
[22, 23]. Thus its mass could be close enough to the D*D thresholds for the resonant coupling 
of the Xci(2P) to the charm mesons to affect the line shape in this region. Its effects on 
the line shapes within about an MeV of the threshold would however be negligible. To be 
more precise, these effects are accurately taken into account through the value of the inverse 
scattering length 7 rc + «7i m - This follows from the universal behavior of an S-wave threshold 
resonance which makes it insensitive to the mechanism for the resonance. If the mass of 
the x c i(2P) is extremely close to the threshold, it is transformed by its resonant couplings 
to the charm mesons into a charm meson molecule whose constituents have a large mean 
separation. Thus far there has been no quantitative analysis of the effects of the xa(2P) on 
the line shape of the X(3872). 

In Ref. [15], Hanhart, Kalashnikova, Kudryavtsev, and Nefediev analyzed the line shapes 
for the X(3872) using a generalization of a parametrization of the line shape for a near- 
threshold resonance proposed by Flatte [24]. Their expression for the line shape can be 
written as |/hkkn(-^)| 2 , where /hkkn(-E') is the scattering amplitude 

/hKKn(E) = (2/ 9l )E f - iT(E)/ 9l + k(E) + (g 2 /g 1 )K 1 (E) - (2/ 9l )E> (3?) 

k(E) is given by Eq. (35a) with r^ = 0, Ki(E) is given by Eq. (35b) with r*! = 0, and 
T(E) is the energy-dependent partial width for short- distance decays of the X(3872). This 
lineshape was also used in a recent analysis by Zhang, Meng and Zheng [17]. Near the 
D*°D° threshold, the scattering amplitude in Eq. (33a) reduces to the universal expression 
in Eq. (2) with 7 re = — [2Ef + <72^(0)]/<7i an d 7i m — F(0)/gi. The coefficient —2/gi of 
the term —(2/gi)E in the denominator of Eq. (37) can be identified as r s /2, where r s is 
the effective range. In Ref. [15], Hanhart et al. found that the best fits to the Belle and 
Babar data are in a scaling region of the parameter space in which the (2/ g\)E term in the 
denominator is negligible. Thus the only relevant parameters are Ef/gi, g 2 /gi, and T(0)/g\. 
The scattering amplitude in Eq. (37) with the (2/gi)E term omitted can be obtained from 
the coupled-channel scattering amplitude f 00(E) in Eq. (33a) by replacing k±(E) in both the 
numerator and the last term in the denominator by /ti(0). Thus the line shape |/hkkn(-E)| 2 
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takes into account some of the effects of the coupling between the neutral and charged 
channels. 

One apparent advantage of the line shape |/hkkn(-^)| 2 that is actually illusory is that it 
is an integrable function of the energy E. The product of branching fractions (BB) F for a 
short-distance channel F of X(3872) in the decay B + — > K + + X can therefore be defined 
in the conventional way by specifying the energy distribution for the final state F to be 

— = T[B + ](BB) F \M KKN (E)\ 2 / J dE'\f RKKN (E')\ 2 . (38) 

This definition of (BB) p appears to be more natural than the prescription for an threshold 
S-wave resonance that we introduced in Eq. (8). However the term in Eq. (37) that makes 
the line shape |/hkkn(-E)| 2 integrable is the last term —{2/gi)E. Because the best fit is 
in a scaling region of the parameter space, varying the parameter gi while holding the 
combinations Ef/gi, g-ij g\, and F(0)/gi fixed has essentially no effect on the line shape very 
near the resonance but it does change the interval of the energy E that gives significant 
contributions to the integral of the line shape. The numerical value of (BB) F is therefore 
determined by the value of g\. Thus the definition for (BB)f in Eq. (37) is in fact an 
arbitrary prescription specified by the parameter g\. 

In their analysis of the line shape of X(3872) in the D°D°7r° decay mode in Ref. [15], 
Hanhart et al. made a serious conceptual error. They assumed that the line shape has the 
form 

^oc|/hkkn(^)| 2 ^ 1/2 W. (39) 

The factor 9(E) emphasizes that the line shape was assumed to be zero below the D*°D° 
threshold. This reflects the incorrect assumption that D°D°7r° events can come only from 
the production of D*°D° or D°D*° above the threshold followed by the decay of D*° or D*° 
and not from the decay of a bound state below the threshold. However the mass of the 
X(3872) is about 7 MeV above the D°D°ir° threshold, so there is plenty of phase space for 
the decay of this bound state into D°D°it . Moreover the X(3872) spends most of its time 
in a configuration in which the charm mesons have large separation, so the D*° or D*° in 
the bound state can decay almost as if they were free particles. The conceptual error in 
Ref. [15] was pointed out in Ref. [18], and an analysis that takes into account the decay of 
the bound state was carried out. 

Zhang, Meng and Zheng have recently carried out an updated analysis [17] of the recent 
data from Babar and Belle using essentially the same line shapes as Hanhart et al. They 
repeated the conceptual error of Ref. [15] by taking the line shape in D*°D° to be given by 
Eq. (39), which does not take into account D°D ti events produced by decays of the bound 
state. They determined the location of the poles in the energy E for the scattering amplitude 
/hkkn(E') in Eq. (37). All of their fits had one pole for which the real and imaginary parts 
of E were less than 1 MeV. This is the pole associated with the S-wave threshold resonance. 
Their fits also had a second pole on a different Riemann sheet of the complex energy E 
whose absolute value was significantly larger than 1 MeV. This pole is an artifact of the 
scattering amplitude in Eq. (37) and has no physical significance. 
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VII. SUMMARY 



We have carried out an analysis of the line shapes of the X(3872) in the J/ipn + n^ and 
D°D°Tr° decay channels using the most recent data from the Babar and Belle Collaborations. 
For the signal, we used the line shapes of an S-wave threshold resonance, which differ 
in several crucial respects from the conventional Breit-Wigner resonance. We took into 
account the experimental resolution in the energy distributions using Gaussian smearing 
functions. In the case of the D D ti channel, we also took into account the assumption in 
the experimental analyses that D°D 7r° events near the D*°D° threshold come from D*°D° 
and D°D*°. 

The parameters for the S-wave threshold resonance are the real and imaginary parts of 
the inverse scattering length 7 re + «7i m and a normalization factor (BB) f that depends on 
the decay channel F. A characteristic feature of an S-wave threshold resonance is that its 
line shapes are not integrable functions of the energy. One consequence is that the product 
(BB)f of the branching fractions for the production of the resonance and its decay into the 
final state F depends on the prescription used to define it. Our prescription for (BB)j/^ 7r + 7T - 
is specified by the analytic expression for the line shape in Eq. (8). Our prescription for 
(BB) D°D°n is specified by the analytic expression for the line shape in Eq. (32). The 
parameters for our fits to the Babar and Belle data in the J/tpn + 7r~ and D D ti decay 
channels are given in Tables I and II. 

Because the line shape of an S-wave threshold resonance is not an integrable function 
of the energy, a prescription is required to define the binding energy and the width of the 
X(3872). Our prescriptions for the binding energy Ex and the width Fx are that the pole 
in the amplitude as a function of the complex energy E are at —Ex — iTx/'i- Given the 
values of 7 rc and 7 im , E x and T x can be calculated using Eqs. (3). An alternative pair of 
variables that can in principle be measured directly are the position E max of the peak in the 
line shape and its full width at half- maximum Ff w h m - Given the values of 7 rc and 7i m , E ma , x 
and rf w hm can be calculated by solving Eqs. (4) and (5). The values of —Ex , Tx, E max , and 
rfwhm for our fits to the Babar and Belle data in the J/ip -k + -k~ and D°D n decay channels 
are listed in Tables I and II. 

We carried out two fits to each of the data sets from the Belle and Babar Collaborations, 
one with 7 im = and one with 7 im as a fitting parameter. The best fits to the smeared 
line shapes in the <//■?/> 7r + 7r~ decay channel are shown in Figs. 1 and 2. The best fits to 
the smeared D*°D° energy distributions are shown in Figs. 6 and 7. The line shapes in 
the J/^it + it~ and D°D°tt decay channels and the D*°D° energy distributions for three of 
the best fits are shown in Figs. 3, 4, and 8. The D°D ir line shape has a peak below the 
D*°D° threshold whose position and width are close to those for the J/ipiT + 7r~ line shape. 
It is this peak that should be identified as the X(3872) resonance. The D*°D° energy 
distribution, which vanishes below the threshold, has a peak above the threshold whose 
width is considerably larger than the width of the X(3872) resonance. Thus measurements 
of the position and width of the peak in the D*°D° invariant mass distribution should not 
be interpreted as measurements of the mass and width of the X(3872). 

In our analyses of the D*°D° energy distributions measured by the Babar and Belle 
Collaborations, we took into account the assumption that D°D°-k events near the D*°D° 
threshold come from D*°D° and D°D*°. Even though the D*°D° energy distributions vanish 
below the D*°D° threshold, our analyses of these distributions gave values for the position 
of the X(3872) resonance that were below the threshold. In our analyses with 7 im as a 
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fitting parameter, the best fits were for 7i m = 0, which is the minimum physical value. In 
contrast, the best fits to the Babar and Belle data on <//■(/> 7r + 7r~ gave 7; m = 15.5 MeV and 
12.0 MeV, respectively. The preference for the value 7i m = in the fit to the D°D ii data 
could be an artifact of the simple model for the D*°D° energy resolution given in Eqs. (30) 
and (31). Because this model is questionable, we regard our analyses of the D°D°7i° data 
as only illustrative. They are no substitute for analyses by the experimental collaborations 
that take all the correlated errors properly into account. In a careful analysis, it would be 
better to use the differential decay distribution in Eq. (32) instead of our analytic expression 
for the line shape in Eq. (22). Finally an analysis of the D°D°n° decay channel similar to the 
original Belle analysis in Ref. [10] would be preferable to one in which D°D°7r° events near 
the D*°D° threshold are interpreted as D*°D° or D°D*°. If such an analysis gave resonance 
parameters for the X(3872) that are close to those from analyses of the J/ipn + n^ channel, 
it would go a long way towards solidifying a consensus in the high energy physics community 
on the nature of the X(3872). 
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